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We consider quantum systems which interact strongly with a rapidly varying environment and 
derive a Schrodinger-like equation which describes the time evolution of the average wave function. 
We show that the corresponding Hamiltonian can be taken to be Hermitian provided all states 
are rotated using an appropriate unitary transformation. The formalism is applied to a variety of 
systems and is compared and contrasted with related results describing stochastic resonances. 
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I. INTRODUCTION 

The study of quantum systems which interact strongly with their environment often presents serious challenges 
due to the possibility that these interactions cannot be neglected and, in addition, also vary rapidly and randomly 
with time [1]. The effects of such external fields is often unavoidable and interesting and can lead to unexpected 
phenomena such as, for example, those studied under the blanket term of stochastic resonances [2]. In this paper we 
will study one subset of such systems. 

We will assume that the interactions with the environment are described by a time-dependent contribution to the 
Hamiltonian denoted by H'{t) which cannot be treated perturbatively. In analogy with a similar problem in mechan- 
ics [3] we will assume self-consistently that the states of the system can be decomposed into a sum of slowly varying 
modes and high-frequency components of small amplitude. Using this decomposition we will show that the time 
evolution of the slow modes is determined by an effective Hamiltonian which, to leading order, depends quadratically 
on the external interaction H' . The formalism assumes that the time scales associated with the interactions with the 
environment are much shorter than all other frequencies in the problem. Denoting by SI a typical frequency of the 
interaction H' , we will obtain a solution to Schrodinger's equation as a series in l/Sl. 

The effective Hamiltonian describes the average time evolution of the system and can exhibit resonances under 
some special circumstances which will be illustrated using simple examples. It is also worth noting that the same 
formalism can be applied to any system evolving according to a Schrodinger-like equation assuming that the operators 
corresponding to the Hamiltonian contain terms which vary rapidly in the evolution parameter. We also provide 
examples of this type of generalization: using geometrical optics we study light-ray propagation in a random media, 
and, we determine the effects of a time-independent potential which varies rapidly with position on the wave functions 
of a non-relativistic particle. 

The paper is organized as follows: in section II we give a description of the formalism and find the effective Hamil- 
^ tonian that will be used in the applications. The behavior of the effective Hamiltonian under unitary transformations 
is studied in section III; the formalism is then applied to various illustrative examples in section IV. In section V we 
give an alternative view of the problem in terms of the Fokker-Planck equation and the results are then compared and 
contrasted with the formalism used in deriving the standard stochastic resonances (section VI). Paring comments 
and conclusions are presented in section VII. Finally, a mathematical detail is relegated to the appendix. 



X 



■ II. QUANTUM SYSTEMS WITH RAPIDLY- VARYING EXTERNAL FIELDS 



We consider a generic quantum system with a Hamiltonian of the form 

H^Ho + H\ (1) 

where H' is time dependent with characteristic frequencies assumed larger than all the other energy scales in the 
system (we take units where h — 1). In general we will also allow Hq to vary with time, but with the restriction 
that the time scale(s) associated with Hq are much smaller than those associated with the time variation of H' . In 
addition we assume that H' is larger than Hq so that, symbolically 
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Hq, Ho/ Ho < H' < H'/H'. (2) 
More specifically we assume that H' admits a Fourier expansion of the type 

H'=Y, H^e-"-\ (3) 

|a;|>n 

where the sum is over a set of frequencies {co} such that the differences also obey |w — u)'\ > ft. In general we will 
allow the Fourier coefficients H^^ to be time-dependent, but, as for Hq, we assume that the corresponding frequencies 
are small compared to fi. Henceforth "slow" will mean "of frequency <C fi". 

In solving the Schrodinger equation for such a system wc will assume that the wave function can be separated into 
a slowly varying piece tp and a rapidly varying (frequency > Q) "ripple" x of small amplitude, 

*=^ + ^- (4) 

slow fast 

We will then match slow and fast terms, noting that even though |x| is small, x can be large. 

We will assume that all quantities can be written as the sum of a slowly varying piece (containing frequencies <C fl) 
and a fast piece (in special cases one or the other may vanish). It then proves convenient to introduce the following 
notation: for any quantity E! 

(S) : slow part of S. (5) 

For example, (*) = i', {H) = H^, {H') = 0. 

To solve Schrodinger's equation for this class of systems we consider a typical term in (3) and define the expansion 
parameter 

rj^-^ (6) 

(alternatively f] ^ \ J dtH'\); the previous restrictions imply that r] < 1. We will now assume that the wave function 
has an expansion in powers of r]: 

* = + Xl + X2 + • • • , Xn ~ (7) 

This expansion is useful for n <C 1/?], beyond this order we typically obtain n— fold products of slowly varying 
quantities which can generate terms with frequency ~ il, and the separation between slowly and rapidly varying 
terms cannot be maintained. These effects are, however, subdominant since their amplitude are suppressed by a small 
factor cx r/" ^ 1. 

Substituting H from (1) and \1/ from (4) in Schrodinger's equation 

= iif, (8) 

we find, to lowest order in 1/fl 

ii) + iXi + 0{nr}^) = Hoi> + //> + 0{HoV, ^rf), (9) 

whence 

iXi = H'yj, hp = Ha^j. (10) 
Since ip is slowly varying, the first equation can be solved to this order in rj by taking ip constant, namely, 

XI =U^; ^=(- [ dtH'\ = -^-e-^'^*. (11) 

To the next order we write 

* = V + Z^V' + X2 + 0(J73), (12) 

and obtain 
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+ iUij + ix2 = Hoi) + HoUi) + H'Ui) + 0{Q.rf^, Horf). (13) 
Note that H'U contains both slow and fast terms. Using the notation (5) we find 

ii)={Ho + {H'U))i), 

iX2 = {[Ho,U] + H'U - {H'U)) V, (14) 

where the second equation can be solved (to the present order in 77) by neglecting the time variation in Hq and tp. 
To this order the average wave function then obeys a Schrodinger-like equation with an effective Hamiltonian [4] 

H^s = Ho + {H'U) . (15) 

It is easy to see that to this order in rj Hgs is Hermitian, however, to order r/^ we find 

ffeff = Ho + {H'U) - {U i[Ho,U] + H'U)) , (16) 

which is not Hermitian: 

ifeff - H^s^ = {H'U +UH' +[U^,Ho]} + ---. 

= idt{U^)+[{U^),Ho]+--- (17) 

and cqnals, to this order the total time derivative of the operator (U"^)- This property corresponds to the fact that 
there is some probabihty "leakage" of order if' from the slowly varying part of the wave function to the rapidly varying 
ripple. This is to be expected since (|xp) = 0{rf') and is non-zero in general. 

The non-Hermiticity of the effective Hamiltonian for the slowly varying modes frequently appears in expansions 
similar to the one considered here [5] ^. This result can be better understood by considering the behavior of the above 
expansion under unitary transformations to which we now turn. 



III. UNITARY TRANSFORMATIONS 

In this section we determine the behavior of the effective Hamiltonian (16) under unitary transformations. We show 
below that the non-Hermitian piece in (16) is modified under such transformations and, in fact, can be completely 
eliminated. 

For the case of a constant transformation, ^ — > with C = it is clear that H — > C^HC and H^s — > C^HegC. 
If the unitary transformation is time-dependent, however, the result is more complicated. We will concentrate on 
transformations of the form 



(18) 



where F is anti-Hermitian, rapidly varying, and of order 77. The Hamiltonian for the transformed states ^ is 



H = e-^He^ - ie-^dte^ 

= H + [H,F] + h[H,F],F]-iF-'- 



F,F 



F,F 



(19) 



A tedious repetition of the procedure outlined in section 11 gives the following expression for the corresponding 
effective Hamiltonian 



^eff = Ho + {H'U) + '-dt {F{F - 2U)) + i {[Ho, [F,U]]) 

+ 1 {[[Ho,U]M) - m'U) - U[Ho, {F - Uf] ) + 0{Hor,', 



(20) 



^In [5] a non-Hermitian term was found already in the first order, the discrepancy between this result and the one obtained 
here is due to different assumptions concerning the time-dependence and magnitude of the various terms in the Hamiltonian, 
leading to different expansion parameters. 
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For a general choice of F this expression is still non-Hermitian. However for the special case 

F = U + 0{rj% 



(21) 



we obtain 

H^s = Ho + ^ {[H',U]) + i {[[Ho,U]M) - {l^H'U) + 0{Horf ,0.7^"), (22) 

which is cxpHcitly Hcrmitian and, in fact, it is identical to the Hermitian part of (16). It is, of course, always possible 

to return to the original frame using = exp{—U H )^'. 

The wave function in the new frame, ^, has an expansion similar to (7) 

* = ^ + Xi + X2 + • • • , Xn-ri", (23) 

where the slowly varying piece -tjj evolves unitarily in time since (22) is Hermitian (at least to order r]^). Using (11) 
we find 

XI = 0. (24) 

The 0{rf) difference between V' and tp quantifies the probability leak into the rapidly varying sector for the original 
frame. The second order term in the wave fimction X2 cannot be determined without a specific choice for the Oirf') 
terms in F. The specific form of the relation between tl) and tl) can also be understood form the expression for the 
non-Hermitian part of the effective Hamiltonian obtained in (17). 

We conjecture that this procedure can be carried order-by-order in r], but since we will not need these higher order 
corrections we will not pursue this further. The expression (22) is the form of the effective Hamiltonian which will be 
used in the following examples. 



IV. EXAMPLES 



The previous results can be applied to a variety of systems. In this section we will consider 5 such examples. Our 
main concern will be to illustrate a wide range of systems that can be studied using the above formalism 



A. AT-level quantum systems 



In this section we consider a quantum system with a finite number of states. This serves, for example, as a model 

for spin or flavor changes in elementary particles; it also describes the basic physics of nuclear magnetic resonance 
and related phenomena. The most general Hamiltonian for these systems can be expanded in terms of the generators 
of SU{N) which we denote by {A"} and which satisfy 

[A«, A*-] = iCabcX". (25) 

We can then write 

Ho = Y.rX\ H^=Y,9Z,>''' (26) 

a a 

where H^^ are the Fourier coefficients of the rapidly- varying Hamiltonian (see (3)). Note that g°L^* = since H' is 
Hermitian. 

Substituting in (22) we find 

-f^eff = y^tp^Ag, 

a 

„c * 6 fd„c*„b n'^ * 

^ =/ 2. -^^'^^'^^ 2^ -^^CabeCeca+ ^ 3a;'(a; + ^^0 ''' '''' ' ^ ^ 

|w|>n |w|>n lo-Uo-'^n ^ ' 
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where the elhpsis denote higher order terms in ij. 

For the particularly simple case of a two-level system with Cabc = "^^abc and Xa = Oa (the usual Pauli matrices) this 
reduces to 



^eff= |f- E ^^^^ + ---| (28) 



It is clear that there is the possibility for the iJ'-induced terms to generate resonances if we allow the to 
vary slowly in time. To see this explicitly we consider (28) taking for simplicity f^=g^ = 0, and /i =constant. 
Substituting we find 

^i=f\ V2 = 0, = f -E' + 0{l/9?), (29) 

where 

^' = 2 E (30) 
|w|>n 

If is allowed to vary slowly with time the effective Hamiltonian will exhibit a resonance when = E' ii\E'\':^ 

This resonance will lead to large transition amplitude provided varies sufficiently slowly: |/^| ^ \J \P\ (which is 
the usual adiabatic resonance condition [6]). A specific example is presented in fig. 1. 

In fig. 1 we also compare the results obtained using the effective Hamiltonian (28) to those obtained solving the 
Schrodinger equation exactly with initial conditions ^{t = Q) = \j}{t = Q) = 1. It can be seen that the solutions 
obtained using iJcff do indeed describe the average behavior of the wave-function provided rj is sufficiently small. We 
conclude that, at least for the case of a 2-level system, the effective Hamiltonian (22) accurately describes the average 
evolution of the system. 

The presence of noise {i.e., H') can then generate unexpected resonances. The condition for these to occur is, 
qualitatively, 

~ (31) 
we will see later that these resonances are related, but not identical, to the well-studied stochastic resonances [2]. 



B. 5-function comb 



A simple model where the above formalism can be applied and which is also exactly solvable is provided by a 2-level 
system with Hamiltonian 

N 

H = ^^hi5{t-ti-nT)- 0<ti<T, (32) 

n i=l 

where the matrices arc Hermitian. This represents a set of N 5 functions which repeats with period T. This type 
of potential is of interest in signal processing [7] and is also similar to the one used in the study the effect of a laser 
beam on a set of charged particles [8] . 
For simplicity we will assume 

^^ = (a* o)' ^^^^ 

in this case Hq = 0, so that H = H', we will also take ft = 2'k/T. 
In order to obtain the effective Hamiltonian we first construct 

-i dt H' = -i'S2hiQ{t-ti), {0<t<T), (34) 
where 9 denotes the step function. U is then the fast part of this quantity. 
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m\\M' (a) (b) 




FIG. 1. Example of a resonant phenomenon induced by the presence of rapidly varying interactions of large amplitude, 
(a): comparison of l^'p obtained by an exact (numerical) integration of Schrodinger's equation (black jagged curve) with \'ip\'^ 
obtained integrating (14) using the effective Hamiltonian in (28, 29) (light superimposed curve). The top light curve is the result 
of integrating Schrodinger's equation when H' — 0; the dotted vertical line denotes the time at which the diagonal elements 
in Heft vanish, (b): Diagonal element of Hcb; for this example we chose E' = 4.78684 and = E'[2/{1 + t'^) + 1/5]. (c,d): 
Non vanishing elements oiH'- the specific expression used was g^ + ig^ = 39.9567e^^^-^^"-°-^''''^°'*'+ 35.6145e*'«''-^''"+° ''^'"'^^* + 
39.1024ei057.84t-o.732253.^ 4239ei208-99*+2 «462. ^ 29.1863ei353.06t+o.7i9083i ^jjig example, rj ~ 0.03 
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u 



Q{t-ti)-l + 



where the slow part of a quantity is obtained by averaging over the period T. 
Using then (33) and substituting in (22) we easily find 



(0 < t < T), 



(35) 



H,^ = -{[H',U]) + - 



1 
T 



(73 + ' 



(36) 



For this system we have ~ 1/T so that i] = max{|Aj|}; it then follows that (36) will be accurate provided |Aj| <C 1. 

This model can be solved exactly by elementary means. Replacing 5{t — U — nT) by a rectangle of height 1/r and 
width T centered t = ti + nT it is easy to see that in the limit r — > 



*(t,+) = e-^^^*(tr), 

where if =ti±5, 5 ^ Q. It follows that 
In the limit where the Aj are small we obtain 



e-^^i*(0). 



(37) 



(38) 



g-iAw ^-ih.N- 



■e = exp 



i^[Ai,A,] + --- 



i>j 

= exp [-iTH^s + ■■■] 



(39) 



which shows that, at least for small Aj, determines the leading contributions to the average time evolution of the 
wave function. 



C. Geometrical Optics example 



The calculations in the previous sections referred to quantum systems, but it clear that any system whose dynamical 
equations can be cast into a Schrodinger-like form can be treated in the same way. In particular for this general case 
there is no need to require if to be a Hermitian operator. 

An example of this situation is provided by the description of light-ray evolution within geometrical optics [9] . For 
small angles the position and direction of light ray within geometrical optics can be described using a two-component 
vector 



(40) 



where h denotes the height with respect to a reference line and a the tilt (assumed to be small) . Any transformation 
of a light ray can be described using a 2 x 2 matrix [9]. In particular 



translation 



1 X 

1 



refraction :( . / ^ i\/d ^ h (41) 
\{ni/n2-l)/R ni/n2 J ^ ' 

where the translation is by a distance x and the refraction is from a medium of refraction index m to another with 

index n2 and i? denotes the radius curvature of the interface. 

We now assume that R and the index of refraction change smoothly though rapidly with distance. We define 

1 dn 
n dx ' 
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so that the general matrix which transports a ray by a distance 6x is 



where primes denote x derivatives and 



M = 



c 



dx 



R 



(43) 



(44) 



The system then corresponds to a two-level quantum system with "time" x and "Hamiltonian" H = iM . A simple 
application of (15) yields 



„ _ . , 1 

-neflf - M g I I 



K^ = {C\\nn)) = +\{{\nnf{l/R)'), 



(45) 



where we assumed {Q') = (;/)= and we kept only the first corrections induced by the rapidly varying terms ^. Note 
that K"^ can be negative and that it vanishes (at least to lowest order) when Rot n are constant. 

The effective operator which determines the translation over a finite distance X is then (assuming for simplicity 
that K is position independent) 



cosh{KX) 
Ksmh{KX) 



{1/K) smh{KX)\ 
cosh(i\:X) ) 



(46) 



(the K-^ < case is obtained by analytic continuation). This matrix is equivalent to a thick symmetric lens with 
radius of curvature R and thickness d such that 



1 
R 



K 
l-n 



taxi\i{KX/2) 



d= —smh{KX), 



(47) 



where n denotes the index of the lens material. 

Thus, within the approximations inherent to geometrical optics, the effects of a region of rapidly varying index 
of refraction and curvature on light rays are equivalent, on average to those of a thick lens of appropriately chosen 
characteristics. For example, using a thick lens with a high index of refraction, n 1 we find 



K\K\ = - 



2n^ 
Rd' 



(48) 



Conversely, a thick lens can be found which completely cancels the effects of A and this can be used to measure the 
fluctuations in the original system (more specifically, those fluctuations which contribute to Q. 

The above expressions suffer corrections form higher-order terms in the expansion in powers of r]. Using (16) ^ we 
find that the next-order term in (45) is 




fc2 



In 7 



In) 



(In n) 



(49) 



where we assumed for simplicity that ((Inn — (lnn))C) and ([Inn — (Inn)]^) are independent of x (in general the 
averaged quantities may still vary slowly with x). The second order correction is negligible provided \K'^\ ^ 

In addition there are deviations form these predications due to the inherent limitations of geometrical optics (for 
example, it is assumed the light rays lie on a plane, diffraction is neglected, etc.). In neglecting them we have tacitly 
assumed that the scale of all fluctuations is large compared to the wavelength and that all reflection and refraction 
angles are small. 



^Tho second expression for K'^ follows from 2v{\iin)/R = (1/Ji)'(lnn)^ - [(lnn)^/ii]' and {[{Innf / R\') = to first order in 
the rapidly varying quantities. 
^As mentioned previously there is no reason to demand Hermiticity in the effective Hamiltonian for this case 
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D. The noisy Jaynes-Cummings model 



In this example we consider a simplified version of an atom interacting with a photon field, as described by the 
Jayncs-Cummings model [10], with the addition of two types of interaction with an external rapidly varying fields. 
We will show that this problem is also well suited for study using the techniques introduced above. We first assume 
that the external field are coupled to the photons, and then directly to the atoms. We then show that both situations 
arc unitarily equivalent. 

The unperturbed Hamiltonian for this model is 

Ho = uoo' a + ^Qo{a'' + aa+) + €(73. (50) 

Denoting by \n; t) and |n; J.) the states with n photons and atomic spin up and down respectively, the eigenstates of 
Ho are [10] 

|n±) = Cn^\n; i) ± sign(K)c„±|n - 1; T), 

Cn± ^ ^ ( 1 ± y- = ) ; K=- — —, (51) 

with energies 

1 



2^0 



\/T+ni^. (52) 



a. Noisy photon field We will now couple the photons to external sources which vary rapidly with time. The 
interaction Hamiltonian is assumed to be 

H' = ^a^ + Ca. (53) 

Substituting (50) and (53) in (22) yields 

Heff = Fo+Wo(|^n-Im(r^)+0(7?3); ^ = i0. (54) 

The difference between H^g and Hq is, in this case, trivial and can be eliminated by a simple change in the overall 

phase of the states. 

Non-trivial terms may arise at higher orders but this would require calculating the (Hermitian version of the) effective 
Hamiltonian to order 77^. Instead of following this miinspiring approach we consider a different way of introducing the 
interaction with the external fields and then show that the corresponding effective operator corresponds to the order 
77^ contribution generated by (53). To this end we first consider a unitary transformation of the form 

S = exp [a^Cit) - aC (t)] ; i( - ivoC + ^ = 0. (55) 

Then the transformed Hamiltonian is 

ifnew = S{Ho + H')S^ + iSS^ 

= Ho- (cT-C + ^+0 - ^Im(C*C) - ^olCl' 

= Ho+K,^, (56) 

which defines H^^^. This shows that the original system is equivalent to one where the external sources are coupled 
directly to the spin of the atoms through C. From its definition it can be seen that ^ is of order rj so that, substituting 
iJ^g^ in (15), gives the effective Hamiltonian up to and including terms of order 77^^ 77^. The effects of this type of 
interaction are considered in the next paragraph. 

b. Noisy spin interaction We now consider 

i/' = Ca++CV_, (57) 

which to lowest order (using (15)) gives 

ffeff = ffo - Im U*^) (T3; C = i^, (58) 
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corresponding to the non-trivial replacement 



e ^ Ceff = e - Im . (59) 

This replacement also describes the leading (average) effect of (53) provided we identify $, = + 0{r])]. If ^ ~ 
then the modification is indeed of order ry^. In order to go back to the original problem we act with S on the states. 

The effects of the external sources is, to lowest order, summarized by the simple shift (59) which corresponds to a 
change in the energy gap between the two spin states of the "atom" of this model. In particular, for noise of sufficiently 
large amplitude, we can have Ces = a;o/2 in which case the Rabi frequency vanishes, En+ = En- and photon number 
is conserved. 

The shift in e can also lead to resonant behavior between states of different n. For example the energies for the 
states |0— ) and |1— ) are equal provided 



1 1 



(60) 



which has a solution only for < wq/^o < 1/2. 

These results are reliable provided ?7 <C 1 which corresponds to |^?| <C 1 and |a;o|, \i^o\, |e| ^ ^- As in the two level 
system resonances occur when \Ho\/Q ^ rj'^ . 

E. Quantum system in an inhomogeneous potential 

The ideas presented in the previous sections can be translated to the case of a particle whose Hamiltonian is of the 
form 

H = -^V^+Vo + Vi, (61) 
2m 

where Vi varies rapidly with position. For this case we consider the time-independent Schrodingcr equation H"^ = E"^ 
and look for solutions 'i' — tp + x where x is of small amplitude but exhibits rapid variation in position while ip is 
slowly varying and of large amplitude. Substituting this Ansatz we find 

.l-V^ + Vo^i^^Ei,-{Vix), 
which can be solved to lowest order giving 

In this case, for any quantity A, (A) denotes the part of A (if any) which varies slowly with position. 

In one dimension the same result can be obtained by converting the time-independent Schrodinger's equation to 
a first order equation for the vector {"if , —id'i/ / dx) and substituting in (22) or (16), and using x as the evolution 
parameter. 

The additional term in Hcn is negative definite and will tend to bind the particle. In particular, taking Vo = 
and assuming Vi vanishes at infinity the effective Hamiltonian H^e will always exhibit a bound state (of zero angular 
momentum) in < 2 dimensions, that is, in < 2 dimensions a rapidly varying potential of zero average will always 
exhibit localized states. The same is true in higher dimensions provided the amplitude of Vi is large enough 

V. PROBABILISTIC CONSIDERATIONS 

In this section we provide an alternative view of the problem using the Fokker-Plank equation. For simplicity we 
consider the case of a two-level system with Hamiltonian 
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H = Ho + H', 

Ho = y^/°o"o, 

a 

H' = J2G''(Ta, (64) 

a 

where {G°} are stochastic variables whose probability function will be described below. 

We will study the Fokker-Plank for the polarization vector tp^atp whose probability density is given by 

V{r,t) = {5^''>{i;\t)am-r))^, (65) 

where the symbol {■■ ■)q denotes the average over the stochastic variables G". In terms of a functional integral we 
will use 

(^)g = / ^ I / / dt'Y,G''{t)Kab{t,t')G\t') \ , (66) 

a V ab ) 

with /C symmetric {JCab{t,t') = JCba{t',t)) and positive definite. We denote by a the inverse kernel /C~^, 

fdsY, ^ac{t, s)a''\s, t') = 5l5{t - t'). (67) 

c 

It is easy to see that a''^{t,t') = {G''{t)G^{t')) . 

We will now restrict further considerations to cases where G°'{t) is correlated with G^{t') only for t close to t' , that 
is for cases where aab{t, t') vanishes except when t ^ t'. In this case we define 

^ab{t) = f dt' aab{t,t'), (68) 

J — OO 

and, following the standard derivation of the Fokker-Plank equation [11,12], we obtain, 

\ a ab / 

where the determine Hq in (64) and La, a= 1, 2, 3 denote the usual angular momentum operators in 3 dimensions. 
There are corrections to this equation but these can be ignored provided a{t, t') is sufficiently localized around t = t'. 
We will now restrict ourselves to situations where a takes the form 

a"'' = ^Dd"'' -\Y1 C^O) 

c 

The first term is commonly used in treating this type of problems, the second term implies a correlation between 
and G^ with a ^ h and is usually assumed to vanish; we will find, however that it is precisely this term that is 
responsible for the resonances described previously. 
Substituting (70) in (69) yields 

= 2 [(f - a) • L - iDL'^] V. (71) 

It is important to note that this choice still corresponds to a positive definite kernel K so that (66) is well defined. 
In order to relate these expressions to the ones obtained previously we first write a in terms of the two-point 

correlator, 

a''^{t,t')= f dt' {G''{t)G^{t'))^. (72) 

J —OO 

We then expand in a Fourier series, 

G"(t) = ^G^ie-"^* (73) 
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(not necessarily restricted to > fl) and assume that (G^G^,)^ ~ for ui + ui' ^ (which is equivalent to assuming 
that the correlator a°'^{t,t') is non-zero for t ^ t' only). In this case 

^ = }l^^'Y.p^^^^''^-^))G- (74) 

UJ 

In obtaining these expressions we included a convergence factor e''*, (5 — > in (73). Note that D will vanish unless 
the G^ are continuously distributed in an interval containing a; = 0. 

Comparing this result to (28) we find that the term f — a in (71) corresponds to the leading term in H^a- In 
addition, however, the Fokker-PIank equation contains the non-Hermitian diffusion term —2iDL'^V which forces V 
to decrease exponentially in time for all but the zero-angular momentum modes [12]. Qualitatively this implies that 
for large times the polarization vector will end up uniformly distributed and the effects of the f — a term will be 
completely washed-out. For intermediate times, however, the presence of the f — a term can lead to interesting effects. 

The general solutions to (71) can be obtained in terms of spherical harmonics ^ Yj^, 

'P = Y. C We-'™^'+''*-''"'"*>l(?), (75) 

Im 

where the exponential is introduced for later convenience. Since V is real the coefficients obey 

4* = (-l)'"d™. (76) 

The above expansion can be substituted into (71) leading, for each I, to a set of 2Z -I- 1 coupled ordinary differential 
equations in t that can in principle be solved for any choice of f, a and D. 

To illustrate this procedure we will consider a special case which is similar to the one often studied when considering 
stochastic resonances. We take 

f - a = i?7oZ + ^Mcos(woi)x, (77) 

corresponding, for example, to — = (so that = = 0), = 0, = ucos{uot) and = —vq. The 
equations for the coefficients in (75) then become 



il + -ucos{ojot) [Vl{l + l)-m{m-l)e''"^'em-i + Vi(/ + 1) - m(m -h l)e-^''°*C+iJ = 0, (78) 

For 1 = the solution is simply =constant which is determines the normalization of V; the equations for Z 7^ 
can be solved numerically using standard techniques. The case / = 1 is of special interest since the coefficients 
determine the average polarization of the system as a function of time: 

(V^(t)<TV'(t)> = J dH{Vv) I J dHv 

= [e-,{t)e'^°'] ,Im [iUt)e'^^'] , ^/,it)^ • (79) 

We will be interested in the possibility of the system resonating at the driving frequency wq- This can be investigated 
by first solving the above equations and then Fourier transforming the result. We define, 

= J dte^-'il{t), (80) 

and we study the behavior of |^J(ci; = wo)] as a function of vq for various values of ujq (note that ^ also depends 
explicitly on vq since this parameter appears in the differential equation (78)). The result is presented in Fig. 2 
which clearly shows an enhancement in the Fourier coefficients of frequency vq when vq = wq, the shape of the curves 
are characteristic of resonant behavior. These resonances are also illustrated by the behavior of (^ip^a^tjj) for various 
values of vq, ojq. An example is plotted in fig. 3. 



4 



It is easy to see that V is, in fact, independent of |r|. 
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<z> 




FIG. 3. Resonant behavior of the z component polarization vector (79), (^ij)'^ a^ij)^ = (i), for a two level system. The 
parameters chosen were u = 0.1, = l/\/3i D = 0.0025, a;o = 5 and, for the case wo / wo, to = wo ± 0.2 
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The same system can be studied using the time averaging procedure of sections II-III provided we assume that 
the corresponding restrictions on the parameters are satisfied. Assuming this is the case, the effective Hamiltonian 
corresponding to the choice (77) is readily seen to be, to lowest order in t], 

( wo/2 cos(woi)\ /SIN 
^•'ff- VcosM -V2 )■ ^^^^ 

The corresponding Schrodinger equation for the slow modes .ffeffV' = ^V' can be solved numerically using standard 
techniques and the solutions arc seen to exhibit resonances whenever vq is an integral multiple of cjq • 

There are, of course, differences between the solutions to the Schrodinger equation associated with (81) and the 
solutions derived from (78). The quantities G° used in obtaining (78) are assumed to be stochastic variables whose 
distribution is determined by (66, 67, 68, 70). In contrast, when deriving (81) we assumed the g'^ are non-zero only 
for > ri, and we also required ^ |wo|, \vo\, \gu>/gu>\- 

The values for a obtained in cases become identical if we assume that the average over the stochastic variables in 
the first case give the same result as the average over time intervals much larger than \/^ in the second case. 

The diffusion term D will vanish, as mentioned above, unless the G^^ are distributed continuously around w = 0. 
This term corresponds to the non-Hermitian contribution (17) which, for this case is simply proportional to the unit 
matrix. Taking /^'"^ = g^ = {], in the example of section IV A (for the case N = 2) we find 

(82) 

which is of order 77^. Note that d = if gl;^ are time independent; this corresponds to the vanishing of D should G^^ 
vanish when oj lies in an interval around a; = 0. 

In concluding this section we note that it is possible to generate the required correlation between g^ and g"^ by 
mixing and filtering two uncorrected functions n^'^. The details are presented in the Appendix. 



eff 



eff 



-id, 



d=4- E 
dt , 



gL + igi 



CO 



VI. COMPARISON WITH THE STANDARD STOCHASTIC RESONANCES 

The resonances described above are reminiscent of the well-studied stochastic resonances [2] that are characterized 
by an increased sensitivity to small perturbations when noise of an optimal amplitude is introduced. This feature 
is also observed when the condition (31) is satisfied, and is illustrated by the behavior of the system studied in the 
previous section. 

More specifically, stochastic resonances occur when there is a match between a noise-induced transition rate, and 
the one produced by an external perturbation. If the latter is assumed to be harmonic of frequency ujo, then typically 
resonances occur when luq ~ tttn- This can be understood by considering a system that initially has 2 degenerate 
minima, such that the harmonic perturbation will first favor one and then the other (alternating with period tt/luq)- 
If the resonance condition on tn is realized, then the times at which one minima is disfavored will coincide with the 
times at which noise-induced transitions to the other minima are most probable, and this enhances the response of 
the system to the external perturbation. This behavior is also observed in the systems studied above, for example, the 
resonances in Fig. 2 occur when vq = too where, according to (81), vq is proportional to the noise- induced transition 
rate. 

There are, however, some technical differences. To illustrate these we consider the following one-dimensional system 
that exhibits stochastic resonances 

X = V'{x, t) + e{t), V{x, t) = Vo{x) + UX cos{uJot), (83) 

where e denotes a stochastic variable, u is a small coupling constant and Vq is a potential with two (degenerate) 
minima. The noise is assumed to obey 

{e{t)e{t')) = FS{t-t'). (84) 

Following the same steps [13] described above it is possible to obtain the Fokker-Plank equation for the probability 
density V{y, t) = {6 {x{t) — y))^ and the corresponding average {x{t))^. The Fourier coefficient of {x{t))^ corresponding 
to frequency loq has an amplitude proportional to {\/F)/ ^Jlo'^ + where A denotes the noise-induced hopping rate 
(the Kramers' rate [14]), InA oc —1/F. For fixed loq this amplitude also displays an enhancement at a certain value 
of F [2] . Comparing these results with those obtained in the previous section we note that 
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• The usual stochastic resonances occur for uncorrclatcd noise obeying (84) while resonance behavior in (64) 

requires the correlations implied by having a 7^ in (70). 

• The resonances described above have the usual shape (see Fig. 2) for the resonant curve. This is not necessarily 
the case for the stochastic resonances usually discussed in the literature. 

• Usual stochastic resonances occur whenever the driving frequency is about half the Kramers' rate, which depends 

exponentially on the noise level F. For the case presented in this paper resonances occur when the driving 
frequency is ~ |a| as defined in (74), and is proportional to the square of the amplitude of the stochastic 
variables G". 



In this work we propose a formalism which makes possible to study systems under the influence of rapidly-varying 
external fields (which are not necessarily perturbative) whose typical frequency we denoted by fl. The formalism 
provides a solution as a power series in and assumes a clear separation of fast (frequencies >ri) and slow 
(frequencies <C fl) modes. 

We showed that the evolution of the slow modes is determined by an effective Hamiltonian which is not neces- 
sarily Hermitian; a point noted in other related calculations [5]. The non-Hermitian contributions to the effective 
Hamiltonian, however, can be eliminated by performing an appropriate unitary transformation. 

The formalism was applied to various classical and quantum systems. In some examples we found that the external 
field can produce a resonant behavior in the system. These resonances are related, but not identical, to the stochastic 
resonances studied in the literature [2] In particular the resonant phenomena studied in this paper occur only when 
the interaction with the environment involves several correlated terms. 

In one particular application of the formalism we argued that the presence of a random time-independent potential 
will necessarily generate bound states in systems of dimension 1 and 2, and in other dimensions as well provided 
the amplitude of the potential is sufficiently large. The connection of this result with the phenomenon of Anderson 
localization [15] are tantalizing and will be considered in a future publication. 



We would like to thank W. Beyermann, R. deCoss and T.J. Weiler for illuminating comments and insights. This 
research was supported in part by US DOE contract number DE-FG03-94ER40837(UCR) and by Conacyt (Mexico). 



In this appendix we describe a simple construction which generates stochastic variables g"" satisfying (70,72), in 
terms of a set of uncorrelated variables n', specifically, we assume 



VII. CONCLUSIONS 
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(85) 



and search for new variables g"" satisfying 



{9%t)g\t')) = -F£{t 



e„bXO(t-f'), 



(86) 



where £ is an even function of its argument while O is odd and satisfy 




(87) 



In terms of the Fourier transformed quantities 




(88) 



15 



we require 



where 



{ni{uo)nj{uj')) =-k5{w + ui')Didij, 



£{uj)* =£{-Lo) = +£ioj), 

d{uj)* = d{-uj) = -d{uj). 



We look for a linear relation between ^" and fii, namely 



Writing 



{Q+ + Q-)i 



with Q±{±oj) = ztQ±(oj) and assuming = ±Q± (where T indicates the transpose) we find 

(?+ + Q- = 1, {Q+,Q-}ij = eikjv'^, 
where u'^ = a'^0/{F£). These equations are solved, for example, by choosing 

Q+ = coshu (1 + a a) , {Q-)ij = sinhueifejfi'^ 



(89) 



(90) 



(91) 



(92) 



(93) 



(94) 



with sinh(2M) — 

It follows that given a set of uncorrclated variables rn it is possible to generate the desired correlated quantities g"' 
through a linear filter defined by the (frequency-dependent) matrix K. 
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